A quantum algorithm for the quantum 
Schur-Weyl transform 

By 

Sonya J. Berg 
B.S. (UC Santa Barbara) 2003 
M.A. (CSU Sacramento) 2005 

DISSERTATION 

Submitted in partial satisfaction of the requirements for the degree of 
DOCTOR OF PHILOSOPHY 
in 

MATHEMATICS 
in the 

OFFICE OF GRADUATE STUDIES 
of the 

UNIVERSITY OF CALIFORNIA 
DAVIS 

Approved: 



Greg Kuperberg (Chair) 



Bruno Nachtergaele 



Jesus De Loera 



Committee in Charge 
2012 



Contents 



1 Introduction 1 

2 Hopf algebra representation theory 6 

2.1 Introduction 6 

2.2 Hopf algebras 7 

2.3 The quantum group of g 1(d) 11 

2.4 The representation theory of Hopf algebras 13 

2.5 Gelfand-Tsetlin type bases 16 

3 The combinatorics of Young tableaux and insertion algorithms 19 

3.1 Introduction 19 

3.2 The combinatorics of Young tableaux 19 

3.3 Insertion algorithms 24 

3.3.1 RSK insertion 24 

3.3.2 Quantum insertion 27 

4 The representation theories of the quantum group U q (d) and the 
Hecke algebra H q {n) 31 
4.1 Introduction 31 



ii 



4.2 The representation theory of U q (d) 32 

4.3 The representation theory of the Hecke algebra H q (n) 34 

4.4 Schur-Weyl duality 36 

5 The Pieri and Schur-Weyl transforms 39 

5.1 Introduction 39 

5.2 The Pieri transform 40 

5.3 The Schur-Weyl transform 43 

5.4 Pieri and Schur-Weyl in the crystal limit 48 

6 A quantum algorithm for the quantum Schur transform 51 

6.1 Introduction 51 

6.2 Quantum probability 52 

6.3 Quantum algorithms 54 

6.4 Time complexity of Schur-Weyl 56 

7 Conclusion 58 
References 60 



iii 



A quantum algorithm for the quantum Schur-Weyl transform 

Abstract 

We construct an efficient quantum algorithm to compute the quantum Schur-Weyl 
transform for any value of the quantum parameter q G [0 , oo] . Our algorithm is a 
g-deformation of the Bacon-Chuang-Harrow algorithm [1], in the sense that it has 
the same structure and is identically equal when q = 1. When q = 0, our algorithm is 
the unitary realization of the Robinson-Schensted-Knuth (or RSK) algorithm, while 
when q = oo it is the dual RSK algorithm together with phase signs. Thus, we 
interpret a well-motivated quantum algorithm as a generalization of a well-known 
classical algorithm. 
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Chapter 1 
Introduction 

This thesis addresses some problems in quantum computation that are motivated by 
quantum algebra. 

The quantum Fourier transform for a finite group G plays a central role in the 
theory of quantum algorithms. This is another name for the Burnside decomposition 
of the group algebra of G, 

V 

which is an isomorphism of Hilbert spaces as well as an isomorphism of algebras. 
Since the Burnside decomposition is a Hilbert space isomorphism and therefore a 
unitary operator, one can ask when it can be expressed by a small quantum circuit, 
or equivalently, when it has a fast quantum algorithm. 

Polynomial-time quantum algorithms for the Burnside decomposition are known 
for many finite groups (see for example [16], [2], [12]). In especially favorable cases, 
the quantum Fourier transform for G yields an algorithm for the hidden subgroup 
problem for G or other groups related to G. In particular, the Shor-Simon-Kitaev 
algorithm (see [16], [17], [9]) to find periods or compute discrete logarithms in any 
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finitely generated abelian group is based on the quantum Fourier transform for finite 
abelian groups. 

The Schur-Weyl decomposition is another transform which is related to the Burn- 
side decomposition for the symmetric group S(n). Given a Hilbert space or qudit 
V = C d , the Schur-Weyl decomposition is 

V® n ^@R x ®V x , (1.0.1) 

Ahn 

where R x is an irreducible representation of the symmetric group S(n), acting by per- 
muting tensor factors, while V x is an irreducible representation of the unitary group 
U(d), which acts simultaneously or diagonally on all of the factors of V. The fact 
that R x is the multiplicity space of V x and vice versa is known as Schur-Weyl duality. 
Recently, Bacon, Chuang, and Harrow presented an efficient quantum algorithm to 
compute a basis refinement of this decomposition [1]. 

In this thesis, we clarify and generalize the Bacon-Chuang-Harrow (or BCH) al- 
gorithm. First, the Schur-Weyl decomposition has a generalization that depends on 
a parameter q from quantum algebra. We replace the unitary group U(d) with the 
quantum group U q (d) = U q (gl(d)), and the symmetric group S(n) with the Hecke 
algebra H q (n). Then the Schur-Weyl decomposition still exists for every q G C which 
is not a root of unity, but the specific linear isomorphism expressed by equation 1.0.1 
depends on q. (If q is a root unity of order r, then the decomposition still exists, 
but it degenerates into a different form when r = 0(n + k).) When q is real and 
positive, then both sides of 1.0.1 are naturally Hilbert spaces and the isomorphism 
is still unitary. 

Our main result is the following theorem which appears in Section 6.4. 
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Theorem 1.0.1. There is an efficient continuous family of quantum algorithms for 
the quantum Schur-Weyl transform for each q G [0, oo]. When q = 1, the algorithm 
is the Bacon- Chuang- Harrow algorithm. The algorithm continuously extends to q = 
and becomes a unitary form of the Robins on- Schensted-Knuth (RSK) algorithm ([15], 
[10]) together with phase signs. The algorithm also continuously extends to q = oo, 
and becomes the dual RSK algorithm without any phases. 

Note the double use of the word "quantum", referring to both quantum com- 
putation and quantum algebra. Those constructions in quantum algebra that are 
non-unitary have no quantum computation interpretation, while many constructions 
in quantum computation only have a pro forma interpretation as quantum algebra. 
Theorem 1.0.1 properly lies in both topics. In fact, the two senses of quantumness are 
slightly incongruous. In quantum algebra, the q = 1 case is called classical or non- 
quantum, because it is the case in which quantum groups become ordinary groups. 
But as an algorithm, the Schur-Weyl transform is not classical when q — 1; it be- 
comes classical when q = instead. The limit q = is called the crystal limit in 
quantum algebra. 

Like the BCH algorithm, our algorithm has running time polynomial in the num- 
ber of qudits n, the size of the qudit d, and loge -1 , where e is the desired accuracy. 
The bound on running time is also uniform in q, assuming that q itself can be com- 
puted quickly. Therefore, our algorithm is efficient in the sense that it is polynomial 
in the number of qudits, for any fixed size of qudit. We do not know whether there is 
an algorithm which is jointly polynomial in n and logrf, i.e., polynomial in the input 
qubit length n(\ogd). 

Our algorithm can be compared to quantum straightening algorithms [11]. Our 
algorithm can be called a Schur-Weyl straightening algorithm, but we emphasize 
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a different interpretation. Straightening algorithms are traditionally interpreted as 
algorithms in symbolic algebra or numerical analysis. As such, the input is not a 
linear number of qubits or qudits, but rather an exponential list of components of a 
vector in a vector space such as V® n . One can make the same distinction between 
a quantum Fourier transform and a classical discrete Fourier transform, which can 
be algebraically the same, but are interpreted differently as computer science. One 
interesting connection between the two interpretations is that a polynomial-time al- 
gorithm for a quantum transform always yields a quasilinear-time algorithm for the 
corresponding numerical transform. (The converse does not hold in general.) 

Finally, in our interpretation and proof of 1.0.1, we will be more precise about the 
basis refinement of 1.0.1. The relevant basis of V x , which can be called the "insertion 
tableau" by extension from the q = and q = oo cases, is its Gelfand-Tsetlin-Jimbo 
(or GTJ) basis. The BCH algorithm and our generalization compute their result in 
this basis essentially by construction — the algorithm is built from a subroutine, the 
Pieri transform (which BCH call the Clebsch-Gordan transform), that stays in this 
basis. We will also prove that the algorithm yields the Young- Yamanouchi-Hoefsmit 
(or YYH) basis of R x , up to sign. (Bacon, Chuang, and Harrow state that it produces 
the Young- Yamanouchi basis without proof.) Finally, specific bases of each V x and 
R x do not quite completely determine a basis of the right side of equation 1.0.1, 
because we could still multiply each summand R x ® V x by a scalar, or in the Hilbert 
space case, by a phase. In this sense, the Schur-Weyl transform is not quite uniquely 
determined. 

This thesis is structured as follows. In Chapter 2 we describe the representation 
theory of Hopf algebras. In particular, we focus on the Hopf algebra of interest in 
this thesis, the quantum group U q (d). In Chapter 2 we also investigate the Gelfand 
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Tsetlin type bases for representations, which have properties desirable for quantum 
computation. In Chapter 3 we describe the necessary combinatorics to discuss the 
representation theory defining the Schur-Weyl transform. We also describe the RSK 
algorithm and a generalization, which we call quantum insertion. In Chapter 4 we 
describe the representation theories of the quantum group U q (d) and the type A Hecke 
algebra H q (n), using the combinatorial language detailed in Chapter 3. We end the 
chapter with the formulation of Schur-Weyl duality, which is central to this thesis. In 
Chapter 5 we define the Pieri and Schur-Weyl transforms with an emphasis on their 
connections with insertion algorithms. Finally, in Chapter 6 we give an introductory 
backgroung to quantum probability and algorithms, and present our main theorem 
1.0.1. 
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Chapter 2 

Hopf algebra representation theory 

2.1 Introduction 

Finite groups and semisimple Lie algebras are familiar examples of algebraic struc- 
tures with nice representation theories. Hopf algebras have an algebraic structure 
which generalizes that of both finite groups and semisimple Lie algebras, while re- 
taining the key properties of their representations. In this chapter we describe Hopf 
algebras and the basics of their representation theory. In Section 2.2 we define Hopf 
algebras and the property of cocommutativity. In Section 2.3 we see an example of 
a noncocommutative Hopf algebra which will reappear in subsequent chapters. In 
Section 2.4 we define the representation theory of Hopf algebras and state some of 
the key theorems in their study. Finally, in Section 2.5 we describe a basis for algebra 
representations which is of both algebraic and computational interest. 
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2.2 Hopf algebras 

The material in this section can be found in [8]. An algebra A over a field F has an 
associative multiplication m : A <g> A — > A. It also has a two-sided unit, which can 
be expressed as a map i: F — > A such that ■ a = a = a ■ for all a in A In 
pictures, the following two diagrams should commute: 
(Associativity Axiom) 



A®A®A^A®A 



idcgim 



A® A' 



A 



(Unit Axiom) 




A coalgebra is obtained by reversing all the arrows. Thus we have a coassociative 
comultiplication A : A — > A® A and a two-sided counit e: A — > F, where the following 
two diagrams should commute: 

(Coassociativity Axiom) 



A®A®A^-A®A 



id®A 

A <g> A*r 



A 



(Counit Axiom) 
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If an algebra A also has a coalgebra structure, so that the maps A, e are alge- 
bra homomorphisms and m, i are coalgebra homomorphisms, then it is a bialgebra. 
Then, a Hopf algebra is a bialgebra with a map called the antipode. The antipode 
is a bialgebra endomorphism S: A — > A where the following three compositions are 
identical: 

A ^A®A^tA®A^A 
A^+A®A^A®A^A 
A^^Y^^A 

We will generally consider algebras over the complex numbers C and the real 
numbers R. 

Example 2.2.1. Given a group G, we can form its group algebra C[G] with basis 
indexed by elements g G G. C[G] is in fact a Hopf algebra with coproduct, counit, 
and antipode map defined by 

A(flO = g®g 
<9) = 1 
S(g) = g- 1 

Example 2.2.2. Let g be a Lie algebra over C. Then its universal enveloping algebra 
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U(g) is a Hopf algebra with coproduct, counit, and antipode map defined by 
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A(a) = a®l + l®a 
e(a) = 
.5(a) = —a 

If A is an algebra over C, then it is a *-algebra if it has a map * : A — > A with 
the following properties: 

(a + b)* = a* + b* 
(ab)* = b*a* 

for a,b G A and A G C. If A is a *-algebra and Hopf algebra so that A(x*) = A(x)*, 
then we call A a Hopf *- algebra. 

If is an algebra over R, then Ac = A^ £g>K C is an algebra over C. On the 
other hand a complex algebra Ac may have more than one decomplexification A®., 
even though there is always an obvious algebra inclusion A^ C A. 

Specifying a decomplexification of A is equivalent to choosing a "bar structure" 
a i — y H that satisfies the axioms: 

a + b = a + b Xa = Xa 

ab = ab a — a. 

This is almost the same as a *-structure, the difference being that a bar structure 
does not reverse multiplication. Given a bar structure on A, the real subalgebra A^ 
is the set of self-conjugate elements a = a. Also, if A has both a *-structure and a 



(Xa)* = Xa* 
a = a, 



2.2. Hopf algebras 

bar structure, then we require that they commute, or 
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(a)* = (a*). 

In some cases, such as for C[G], the antipode map is involutory, and the *-map 
is essentially a conjugate-linear version of the antipode map. In other cases, the 
antipode map will not be involutory and there is some other *-map making the 
algebra into a Hopf *-algebra. 

We end by describing the condition of commutativity and define the analogue for 
the coalgebraic structure of a Hopf algebra. Let A be an algebra with multiplication 
map m. Define the flip map r: A® A — > A® Ahy r(x ®y) = y ® x. One way of 
defining A to be commutative is by requiring the following diagram commute: 

A ®A > A ® A 




A 



Thus, if A is a coalgebra with comultiplication map A, we define cocommutativity 
by requiring the following diagram commute instead: 

A® A^ A® A 




A 



Note that both C[G] and U(g) are generally noncommutative but always cocom- 
mutative. In the next section we'll examine a Hopf algebra which is both noncom- 
mutative and noncocommutative. 
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2.3 The quantum group of Qi(d) 

The Hopf algebras C[G] and U (g) we saw in section 2.2 are generally noncommutative 
but always cocommutative. In this section, we introduce an example of a Hopf algebra 
which is noncommutative and nonco commutative: the quantum group. Quantum 
groups as defined independently by Drinfeld [5] and Jimbo [7] are deformations of 
U (g) for g a Lie aglebra. 

In this section we consider the Lie algebra gl(d) which is isomorphic to End(d), the 
set of linear maps on C d . The generators of U(gi(d)) are and fa for 1 < i < d — 1, 
and hi for 1 < i < d. The relations on the generators are called Serre relations and 
are given by: 



[hi,hj] 


= for j ^ i 




= [hi,ej] = 








= [fijj]=0 for ]i — j| > 1 




= ^ (e?e i± i + e i±l ef) 


fifi±lfi 


— 2 {fifi±l + fi±lfi) ■ 



The associated Drinfeld- Jimbo quantum deformation of U(gl(d)) is called a quan- 
tum group, and is written U q (gl(d)) which we will abbreviate to U q (d). The parameter 
q is a complex number not equal to zero or one. The generators of U q {d) are and 
fi for 1 < i < d — 1, and q ±hi / 2 for 1 < z < d. The generators q ±hi / 2 can be inter- 
preted as formal exponentials rather than actual powers of q. The formal notation 
is meant to imply that these generators commute with each other and that q~ hi l 2 is 
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the reciprocal of q hi ^ 2 , using addition in the exponent. In these formal exponentials, 
we also let ki = hi — h i+ \. 

We use the notation [n] for the quantum integer defined by the formula 



\n\ = ^ Q ! = q^ 1 + q n ~ 3 + ■■■ + q'^ + 



q-q 



-i 



Extending the notation to operators, we write [hi] 



q-q- 1 



The relations on the U q (d) generators are g-deformations of the U(g) Serre rela- 
tions, and are given by 



[q h * /2 ,q h i /2 ] = ,for % ^ j 

q x l 2 e^q hi l 2 for i — j 

q~ l / 2 ejq hl / 2 for i — j + 1 

A/ 2 otherwise 



q h ^ 2 e 3 = i 



[€-ii &j 



e j1 

q~ l l 2 fjq hi l 2 for i = j 
q x l 2 fjq h ^ 2 for i = j + 1 
otherwise 
eijj] = 5ij[hi] 

[fi,fj] = Q, \i~j\>2 



[2] 



e 2 i e i± i + e i±1 e 2 



fifi±lfi — "pj (fifi±l + fi±lfi 



Interestingly, the Hopf algebra structure on U(g) can also be deformed so that 
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U g (d) is a Hopf algebra. For example, the coproduct map becomes 
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A(q h ^ 2 ) = q h ^ 2 ® q h ^ 2 

A(ei) = ej ® g" fcl/2 + q k ' /2 <g> e* 

A(/ l ) = /,®g- fcl/2 + ^ /2 ®/,. 

There are other deformations of the Hopf algebra structure that result in different 
coproduct maps. For example, we could replace A as defined above by r o A, which 
is distinct from A by noncocommutativity. 

When q is real and positive, U q {d) has a *-map defined by 

e* =ft ft = e l (q hi/2 Y = q hi/2 . 

This *-map makes U q (d) into a Hopf *-algebra. When q is real and positive, U q (d) 
also has a bar structure in which all of the generators are real, and they generate a 
real Hopf *-algebra U q {d)^_. 

In future sections we'll restrict to the case when q is real and positive so that we 
can use the associated * and bar structures. 

2.4 The representation theory of Hopf algebras 

A representation of an algebra A is a vector space V and a linear map p : A -> End(V) 
which preserves the multiplication and unity, i.e. p(ab) = p(a)p(b) and p(l) = 1. The 
action p can be implied so that p(a)v is written av, or in quantum notation as a\v). In 
the rest of this section, we fix the assumptions that our algebra A is a Hopf algebra, 
and our representations V are defined over C and are finite dimensional. 
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Two representations V and W of A are isomorphic if there exists a linear bijection 
T : V W that commutes with the action of A, i.e., T(av) = aT(v) for all a G A, 
v G V. 

The representation V is irreducible if it has no non-trivial subspaces that are closed 
under the action of A. In this thesis we use the abbreviation irrep. For example, the 
counit map e: A — > C defines a trivial representation, which is irreducible since it's 
one dimensional. 

Given two representations V and W of A, there is a well-defined representation 
structure on the direct sum V © W given by 

a(v © w) = av © aw. 

A representation V is semisimple if it is isomorphic to a direct sum of irreps. 
(Likewise, an algebra A is called semisimple if all of its representations are semisim- 
ple.) The number of occurrences of an irrep W in V is called the multiplicity of W . 
If the multiplicities are all or 1, then V is called multiplicity-free. 

Assume for the moment that A and its subalgebras are semisimple. Given a 
representation V of A and a subalgebra BOA, the restriction of V to B will be 
denoted by Res s y. When V is an irrep of A, Res^V is typically not an irrep of B, 
but by semisimplicity Res^y decomposes as a direct sum of irreps of B. A rule for 
describing the decomposition of Res^l/ into irreps is called a branching rule. If for 
all irreps V of A, the branching rule for ResgV is multiplicity-free, then the inclusion 
B C A is called a G elf and pair. 

Given two representations V and W of A, the coproduct map A : A — > A © A is 
used to define a representation structure on V © W. 

The antipode map is used to define a dual representation. Given a representation 
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V of A, define V* to be the dual space of linear functionals on V. Then, the action 
of A on V* is defined by a(v\ = (v\S(a). 

If A is a Hopf *-algebra, then V is a *-representation if p(a*) = p{a)* where the * 
on the right side is the Hermitian adjoint. (The Hermitian adjoint makes the algebra 
End(y) into a *-algebra.) This generalizes the notion of a unitary representation of 
a group. In particular, a *-representation V is automatically semisimple: If W is a 
subrepresentation of V, then so is its orthogonal complement W^. 

Although such a V might possibly have non-orthogonal irreducible decomposi- 
tions, it always has an orthogonal irreducible decomposition. If V is multiplicity-free, 
then its irreducible decomposition is unique and therefore orthogonal. 

Our analysis so far carries over verbatim to representations of algebras over R. 
Quantum computation is defined over C, and we will ultimately be interested in 
connecting representations over K with representations over C. 

If Vr is a representation of then 

Vc = Mr ®m C 

is a bar representation of Ac- But note that even if is irreducible, Vz may or may 
not be irreducible. If End(Mg) = R, then Vc is irreducible, while if End(Vi) = C 
or End(Vft) = H (the quaternions), then Vc has two irreducible summands. In the 
former case, we will say that Vr is strongly irreducible. 

Example 2.4.1. C[G] is a Hopf *-algebra with *-map defined by g* = g~ l . Note 
that in this case a representation being a *-representation is the same thing as it 
being unitary as a representation of G. Also C[G] has a standard bar structure with 
g = g, so that its decomplexification is the real group algebra M[G]. 
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When C[G] is finite dimensional it has additional properties for its irreps. For 
example, there are finitely many distinct irreps of C[G], indexed by the conjugacy 
classes of G. And we always have semisimplicity of representations of C[G]. 

Example 2.4.2. Every continuous representation V of a connected Lie group G is 
also a representation of the universal enveloping algebra U(g) and it has the same 
subrepresent at ions . 

If 0m is a real Lie algebra and gc is its complexification, then U(qc) has both 
a natural bar structure — where the real subalgebra is C/(jJb) — and a natural *- 
structure. Since U(qc) is generated as a complex algebra by g^, we define these 
structures by letting 

a = a a* = —a 

for a e K . 

2.5 Gelfand-Tsetlin type bases 

In this section, we describe bases for irreps with special algebraic and computational 
properties. We will be interested in a tower of algebras 

C = A cA 1 cA 2 c...cA n 

and we will use the abbreviation 

RfflJ_ 1 V = Re B £_ 1 V 

for the restriction of a representation V of A}.. 

Suppose that each inclusion A^_\ C A n is a Gelfand pair. Then if V = V n is an 
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irrep of A n , Res^^V is a direct sum of irreps V n -\ of A n -i, and by induction each 
Res^_ 1 Vfc is a direct sum of irreps Vk-i of A^-i- As a result, V is expressed as a direct 
sum of irreps Vo of A = C, and all such irreps are isomorphic and 1-dimensional. 
Thus V has a basis of lines which are encoded by flags 

C = V CV 1 C---CV n = V. 

This line basis is called a Gelfand Tsetlin type (GTT) basis. By extension, any vector 
basis that refines the GTT line basis is also called a GTT basis. Note that in the 
encoding, the number of bits a GTT basis vector requires is the sum of the bits 
required to encode each summand 

To get a sense of the significance of a GTT basis, note that whenever a G Ak and 
v G T4, then av G also. This means that we can express the action of an element 
a on V in the setting of a lower-dimensional algebra, which naturally gives rise to 
a recursive structure. However, note that a GTT vector basis of an irrep V is not 
unique; only the corresponding line basis is unique. The computational strength of 
a GTT basis can still depend on how its vectors are scaled. 

Remark 2.5.1. In some articles in quantum computation, if V is an irrep of a group G 
and H C G is a subgroup, then a basis that refines a decomposition of Res^r V is called 
subgroup- adapted. The analogous notion for us is a basis that is subalgebra- adapted. 
In this terminology, a GTT basis is recursively adapted to a tower of subgroups or 
subalgebras. 

If each algebra A^ is a *-algebra and V is a ^representation of A = A n , then a GTT 
basis is automatically orthogonal, because each restriction Res^_ 1 V/ c has an orthogo- 
nal decomposition. We further require that a GTT vector basis of a "^representation 
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be orthonormal, so that the basis is usable in quantum computation. However, even 
when GTT basis vectors are orthonormal, their phases are still not determined by 
the GTT property. 

If V and W are two irreps of an algebra A, with given GTT bases, then technically 
their combinations V © W and V ® W do not have GTT bases. However, we can 
define standard bases by taking the direct sum and tensor bases, respectively. These 
combinations are GTT bases with respect to the action of A <g> A instead. 
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Chapter 3 

The combinatorics of Young 
tableaux and insertion algorithms 

3.1 Introduction 

The representation theories of the algebras described in this thesis are indexed by 
combinatorial objects called Young tableaux. In this chapter we describe the com- 
binatorics of these objects. In Section 3.2 we define Young tableaux and state some 
of their key properties. In Section 3.3 we describe insertion algorithms for operating 
on Young tableau, which will connect to some interesting representation theory in 
subsequent chapters. 

3.2 The combinatorics of Young tableaux 

A partition A is a list of non-negative integers 



A — (Ai, A 2 , . . . , Xd) 
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such that 

Ai > A 2 > ■ ■ • > X d - 

We say that A is a partition of n, or A h n, if J2k Afc = n. The length of A, denoted 
£(A), is the number of non-zero entries of A. 

A partition A has an associated Young diagram, which is a horizontal histogram 
with £(\) rows; the kth row has A^ boxes. 

Example 3.2.1. The Young diagram of A = (3, 2, 1, 1) is 



If n and A are partitions so that the Young diagram of \i is contained in the Young 
diagram of A, then we write /iQ. If A and \i differ by a single box then A is said to 
cover /i. 

When // C A we can form a Young diagram of skew shape given by A \ \i which 
means removing the boxes in the Young diagram of \i from the boxes in the Young 
diagram of A. 

Example 3.2.2. If A = (3,2, 1, 1), and = (1, 1), then the skew shape A\/x is given 
by 



If A \ /i has at most one box in each of its columns, then it is called a horizontal 
strip. 

A Young tableau of shape A (including skew shapes) is a filling of the boxes of the 
Young diagram of shape A with positive integers. If t is a Young tableau of shape A, 
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we write sh(t) = A. We will use special types of tableaux called semi-standard and 
standard Young tableaux. 

A Young tableau of shape A is semi-standard (abbreviated SSYT) if its entries 
weakly increase from left to right and strictly increase from top to bottom. 

Example 3.2.3. An example of an SSYT with shape (3, 2, 1, 1) is given by 



An example of an SSYT with skew-shape (3, 2, 1, 1) \ (1, 1) is given by 



ii 



We denote the set of SSYT of shape A with entries in {1, . . . , d} by SSYT(A, d). 
When the value of d is obvious, we suppress it and write SSYT (A). 

A Young tableau of shape A h n is standard (abbreviated SYT) if its entries are 
in {1, . . . , n} and strictly increase both from top to bottom and from left to right. 

Example 3.2.4. An example of an SYT with shape (3,2, 1, 1) is given by 



1 


2 


4 


3 


5 




6 






7 







An example of an SYT with skew-shape (3, 2, 1, 1) \ (1, 1) is given by 
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We denote the set of standard Young tableaux of shape A by SYT(A). 

If v is a horizontal strip, then the SYT u obtained by filling the Young diagram 
of v with letters from left-to-right is called ordered. (This is not in general defined 
for skew-tableau, but exists for horizontal strips.) 

Example 3.2.5. The ordered SYT of the horizontal strip v — (3, 1) \ (1) is given by 

^13" 
U = \l\- 

Given any t e SSYT(A, d), let be the restricted tableau in SSYT(A, k) obtained 
by removing all boxes from t with numbers larger than k. 

Example 3.2.6. Let 



Then 





t -- 
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1 


2 












3 












1 


1 


2 






) = 


1 


1 



Note that the skew shapes sh(t ( - k ^)\ sh(t^ k ~^) for an SSYT t are always horizontal 
strips for each k, so we give them the label A^(t) where sh(t) = A. 

Finally, the residue of a box b in a Young tableau t is the difference of its coor- 
dinates. In other words, if b has coordinates then its residue is res(b) = j — i. 
Two boxes in a Young tableau have the same residue if and only if they lie on the 
same diagonal. 

The axial distance a between two boxes b and b' is defined to be the difference in 
their residues, given by a = res(b) — res(b'). It's described as a distance because it 
counts the number of boxes in any path in the Young diagram from the box b to the 
box V where moves left and down count for +1 and moves right and up count for — 1. 



3.2. The combinatorics of Young tableaux 



23 



Example 3.2.7. Consider the SYT 



1 



2 3 



t = 



4 



5 



6 



The axial distance from the box containing three to the box containing four is 3. 
Note the axial distance is antisymmetric, so the distance from the box containing 
four to the box containing three is —3. 

Given a horizontal strip A, we define for i < j to be the axial distance from 
the last box in the zth row to the last box in the jth row. (Note we can equivalently 
use the boxes in the second-to-last, third-to-last, etc., positions, and get the same 
values.) 

We will see the axial distances in subsequent chapters when we describe matrix 
coefficients that derive from representation theory. One main reason these axial 
distances are chosen is that they sum in a very natural way: 



Example 3.2.8. If v is the shape given by 



then a 12 = 3, a 23 = 4, and ai 3 = a 12 + a 23 = 7. 
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3.3 Insertion algorithms 

Given an SSYT t and a letter i, we will add a new box to sh{t) to make space for 
an extra letter and insert the % into the tableau t, possible rearranging other letters 
in the process. In this section we review the well-known RSK insertion algorithm 
and introduce a generalization that we call quantum insertion. The RSK algorithm 
can be found in [10] and [15], while quantum insertion is our way of describing the 
techniques found in [3]. 

3.3.1 RSK insertion 

The first insertion algorithm we examine is called Robins on- Schensted-Knuth (abbre- 
viated RSK), and denoted (i — — >■ t). This insertion algorithm produces a unique 
tableau, given by the rules: 

1. If i is greater than or equal to all the numbers in the first row of t, then add i 
to the end of the first row of t. 

2. Otherwise, pick the leftmost box in the first row containing a number j > i. 
Replace j by i. (This process is referred to as i bumping j.) 

3. Repeat steps (1) and (2) for j starting with the second row. Proceed inductively. 
Example 3.3.1. Start with 



1 


1 


2 


2 


3 




3 






4 







If we choose to insert a letter i > 2, such as i = 4, then it will be added to the end 



3.3. Insertion algorithms 
of the first row: 



25 
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3 
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However, if we choose to insert the letter i — 1, it will bump the two out of the first 
row: 
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4 







The two will then bump the three out of the second row, which will itself get added 
at the end of the third row. Therefore, 
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2 




3 


3 




4 







There is also a dual RSK algorithm, denoted (i RSK > t), which can be thought of 
as the standard RSK algorithm applied to columns instead of rows. Thus, the dual 
RSK algorithm also produces a unique output, given by the rules: 

1. If i is larger than all numbers in the first column of t, add i to the end of the 
first column of t. 

2. Otherwise, pick the topmost box in the first column that contains a number 
j >i. Replace j by %. 

3. Repeat steps (1) and (2) for j starting with the second column. Proceed induc- 
tively. 

Given a word w — W\ . . . w n , we can extend the algorithm by induction to define 



RSK KbK / KSK / / K5R, \ \ \ 

w > = w n > {w n -i > (. . . (w 2 > Wi) ... )) 
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RSK 

w — >- 



w' 



Given the word w' = 221 we still obtain the same output as w, i.e. w 
— On the other hand, w" = 122 results in a different tableau 



RSK 



The above example shows that the RSK map is not invertible. However, note 
there is a way of distinguishing w -^A- and w' by a recording tableau which 
tracks the order new boxes are added in the sequence of insertions, as is done in the 
following example. Then we define rsk(w) = P{w) x Q(w), where P(w) = w -^-A 
and Q{w) is the recording tableau. 



Example 3.3.3. If w = 212, 



RSK lU 



whereas if w' = 221, 



rsk (w 



so that rsk(w) ^ rsk(w'). 

The proof of the following theorem that the RSK map is a bijection can be found 
in [10]. 
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Theorem 3.3.4. Let V] 1 be the set of words in d letters of length n. Then, the RSK 
map rsk(w) = P(w) x Q(w) is a bijection between V£ and the disjoint union 

[] SSYT(X) x SYT(X). 

Ahn 
e(X)<d 

3.3.2 Quantum insertion 

In this subsection we consider a generalization of RSK we call quantum insertion, 
or q-insertion, and denoted (i t). Given an SSYT t and a letter i, q- insertion 
produces a set of output tableaux, one of which is (i t). The rules for constructing 
the output tableaux are given by: 

1. In all possible ways add a new box to sh(t). 

2. In all possible ways, take either of the following two steps. 

• Insert i into the new box. If this step is taken, the algorithm terminates. 

• For any letter j > i, i can replace (or bump) j. In this case step 2 is 
repeated inductively with j. 

Example 3.3.5. Let 
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t = 
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3 













and suppose we wish to insert a 2 after adding a box to the second row. 
Then, starting with 
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1 


2 


2 


3 
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we insert a 2. It replaces the 3 because it can't take over the new box, and it can't 
replace the 4. 



1 


1 


2 


2 


2 









We then repeat the procedure with the 3, which can either take over the new box or 
replace the 4. 



2^> 
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2 




4 






3 





We define a bumping sign for an output tableau as follows. For each letter involved 
in the bumping procedure, multiply the bumping sign by a —1 if the letter moves to 
a lower row in the tableau. Note that for RSK, the bumping sign can be ±1 whereas 
for dual RSK the bumping sign is always +1. 

Analogous to the RSK map, given a word w = W\ . . . w n , we extend the q-insertion 
algorithm by induction to define 



w 



3 qlNS / qlNS / / qlNS \ \N 

-> = w n — y (w n _i — y (. . . [w 2 — y wx) ...)). 



Unlike the RSK algorithm, if w' is a permutation of w, then the sets w ^— y and 
w' ^— y are equal. Also unlike RSK, the output of w is an entire set of SSYT, 
and sometimes there is more than one insertion path in w which produces an 
SSYT t. 

As we see in the following example, we can distinguish outputs by attaching a 
recording tableau which tracks the order in which new boxes are added during the 
insertion process. Then we define q ms(w) = {P q (w) x Q q (w)}, where P q {w) is a SSYT 
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in w and Q q {w) is the associated recording tableau. 
Example 3.3.6. Letting w = 122, 
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qlNS 

w >- 



12 2 
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2 





So, the second and third tableaux in w qINS > are equal. However, we can distinguish 
the output tableaux by attaching a recording tableau to each: 



qINS(W) 
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3 



X 



X 



Thus far, the reason for using the word "quantum" in the context of a combina- 
torial insertion algorithm is unclear. In the rest of this section we describe the reason 
for this choice. Much of the material can be found in [3]. 

Define the weighted q-insertion map by 

q iNs(u>) = ^c P , Q P q (w) x Q q (w), 



for a choice of nonzero constants Cp t Q G C[q, q^ 1 ]. 

The choice of coefficients Cp t Q that interest us is determined by representation 
theory and will be described in chapter 5. The connection between RSK and q- 
insertion becomes clear in the following theorem, which can also be found in [3]. 

Theorem 3.3.7. Let Vd be the vector space over C[q, q^ 1 ] with basis {1, . . . , d}, and 

consider the vector space V^f" 1 of words of length n. Let V x and R x be the vector 
spaces over C[q, q^ 1 } with bases SSYT(X) and SYT(X), respectively. 

Then, there exists a choice of coefficients cp^ so that the weighted q-insertion 
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map q iNs(w) = J2 c P,QPq( w ) x Qg( w ) defines a vector space isomorph' 



Xhn 
l(\)<d 
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Chapter 4 

The representation theories of the 



4.1 Introduction 

In this chapter we describe the representation theories of our quantum algebras of 
interest. In Section 4.2 we present the respresentation theory of the quantum group 
seen in Section 2.3. In Section 4.3 we present the representation theory of the Hecke 
algebra H q (n), which is a (/-deformation of the symmetric group algebra C[jf?(n)]. 
These Hecke algebras are not themselves Hopf algebras, but for almost all choices of 
q their representations are isomorphic to those of CfS^n)], which is a Hopf algebra. 
Finally in Section 4.4 we describe the correspondence known as Schur-Weyl duality 
between representations of the quantum group and Hecke algebra. 



quantum group 




and the 



Hecke algebra H q {n) 
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4.2 The representation theory of U q (d) 

We defined the quantum group U q (d) in Section 2.3 as an interesting example of a 
noncommutative and nonco commutative Hopf algebra. In this section we describe 
its representation theory. 

Recall that we restrict the values of q to real and positive in order to make use of 
the star and bar structures available in this case. For these values of q, the irreducible 
representations of U q (d) are isomorphic to those of the unitary group U(d). (This 
is true for other values of q as well, namely those values of q which are not roots of 
unity or zero.) The irreps of U(d), and hence the irreps of U q (d), are in bijection with 
partitions A whose length is bounded by d, written £(X) < d. For the representation 
indexed by A we write V x . 

Restricting to all the generators except e^-i, fa-i, and q hd , we realize a copy of 
U q (d — 1) inside U q (d). The branching rule associated to this pairing is given in the 
following theorem. 

Theorem 4.2.1. The algebras U q (d) and U q (d—1) form a Gelf and pair. In particular, 
ifV x is an irrep ofU g (d), then 

ResUV x = V". 

X\/i horizontal strip 

e(ji)<d-i 

The branching rule (4.2.1) implies that a GTT basis for the irrep V x can be 
written \v t ) with t e SSYT(X). Then V x is the span of the elements \v t ) so that 
(v t | v 8 ) = 5 t)S . 

The specific GTT basis we use is called Gelfand-Tsetlin-Jimbo (GTJ). The for- 
mulas described in the rest of this section can be found in [3]. The action of the 
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generator q hi on the GTJ basis is the easiest to describe and is given in the following 
theorem. 

Theorem 4.2.2. Let \v t ) be a GTJ basis element indexed by the SSYTt. Then, 

q h ^\v t ) = q Xi{t)/2 \vt), 

where Xi{t) counts the number ofi'sint. 

The generator fi acts on \v t ) by turning an instance of % in the tableau t into an 
i + 1 (in all possible ways, i.e. in superposition). In other words, letting \v tk ) be the 
vector indexed by tableau tj, where the last % in row k is changed into ani + 1 but is 
otherwise identical to t, or zero if this is not possible, then 

/iK) = UK I fi I Vt)M 
k 

for some choice of coefficients (v tk \ fi \ v t ), which we call the GTJ statistic. The 
action of on \v t ) is also defined with GTJ statistics using the relation e* = /j. The 
GTJ coefficients are complicated notationally, but in principle derive from simple 
combinatorial properties, in particular axial distances, of the SSYT t. Recall we 
defined the axial distance a^- of t, and associated horizontal strips AW(t) in section 
3.2. 

Theorem 4.2.3. Let \v t ) be a GTT basis element of V x indexed by the SSYT t. 
Then the GTJ statistic is given by 

k i /. i v.) = .k"nAr» + 1] iff hi_f] >^ + Mg+ a (4.2.D 

\\ j= i i a jk\ [a>jk + J- J 
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4.3 The representation theory of the Hecke alge- 
bra H q (n) 

The Hecke algebra H q {n) is a certain (/-deformation of the group algebra C[S'(n)]. 
(More precisely, we consider a Iwahori- Hecke algebra of type A. There are also other 
kinds of Hecke algebras.) Note that the Hecke algebra H q (n) for q ^ 1 is not a Hopf 
algebra, so the results in this section are proved independently of the theorems for 
Hopf algebra representation theory. 

The Hecke algebra H q (n) with complex parameter q has generators {Ti, . . . , T n _!} 
with relations 

TiTj = TjTi for \i - j\ > 1 

Ti rri rri rri rri rji 

iJ-i+l 1 ! — 1 i+l 1 i 1 i+1 

(T l -q- 1 )(T l + q) = 0. 



The first two relations are known as the braid relations and the third is the quadratic 
relation. (We use the generators used by Jimbo [7]; the generators due to Iwahori 
are slightly different.) When 5=1, the third relation simplifies to Tf = 1, so that 
in this case Tj represents the transposition Sj = (i,i + 1) in the symmetric group. In 
other words, Hi{n) = C[S(n)}. 

As with the quantum group U q (d), we restrict to the case when q is real and 
positive. In this case, H q (n) has both a *-structure and a bar structure, defined by 

T* = Ti Ti = 



For these values of q, the irreducible representations of H q (n) are isomorphic to 
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those of C[S'(n)]. (This is true for other values of q as well, namely those values of q 
which are not roots of unity or zero.) The irreps of C[S'(n)], and hence the irreps of 
H q (n) are in bijection with the conjugacy classes of S(n), so are indexed by partitions 
of n. For the representation indexed by A h n we write R x . Importantly, for these 
values of q, H q (n) representations remain semisimple. 

It is known that the dimension of R x equals the number of standard Young 
tableaux of shape A (given by, for example, the hook length formula). Therefore, 
there is a basis of R x indexed by SYT(A). We describe below how GTT bases are 
naturally described by SYT(A). 

The Hecke algebra H q (n) contains many copies of H q {n — 1); we consider the 
one obtained by restricting to the generators T\, . . . ,T„_2- Thus, we can describe 
the restriction of R x to H q (n). The corresponding branching rule is multiplicity-free 
and has a nice combinatorial description in terms of the covering relation of Young 
diagrams. 

Theorem 4.3.1. The algebras H q {n) and H q {n — 1) form a Gelfand pair. In partic- 
ular, if R x is an irrep of H q (n), then 

Res^_ 1 R x = R". (4.3.1) 

A covers /i 

Theorem 4.3.1 implies a GTT line basis with elements indexed by sequences of 
partitions pairwise differing by a single box, i.e. standard Young tableaux. The 
vector basis of R x we use in this thesis that is a refinement of the GTT line basis 
defined by Theorem 4.3.1 we call the Young- Yamanouchi-Hoefsmit (YYH) basis. We 
write YYH basis elements as \r t ) where t e SYT(A). Then R x is the span of the 
elements \r t ) with (r t | r s ) = 5 tiS . 
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Define the following action of S(n) on the basis element r t : 

• If i and i + 1 are in the same row or column of t then r s .. t = 0. 

• Otherwise, r s .. t = r t > where t' is the standard tableau obtained by switching i 
and i + 1 in t. 

The action of H q (n) on the YYH basis defined below in Formula 4.3.2 is a normalized 
version of that given in [14], building on that found in [6]. 

Theorem 4.3.2. Let a be the axial distance in t from the box containing i to the box 
containing i + The action of H q (n) on R x with the YYH basis is defined by 

Ti\r t ) = fHr t ) + \A-OlK-*> (4-3-2) 

[a] v m 

Example 4.3.3. Consider f?( 2,1 ) with basis 

- - ( » : ) 



Then, 




4.4 Schur-Weyl duality 

Let V be any finite-dimensional vector space over C. Then V® n is a representation 
of C[5'(n)] via the simple permutation action 



7r(vi ...V n )= ® ■ ■ ■ ® %-!(„)■ 



(4.4.1) 
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The vector space V® n is also a representation of the Hecke algebra H q (n) via a q- 
deformation of the permutation action defined in 4.4.1. In particular, the generator 
on a vector in V® n by the identity on all factors except the ith and % + 1st 
ones. On these two factors, it acts by 



T\vi)\vj) = < 



\vj)\vi) iii<j, 
(q- 1 - q)\vi)\vj) + if i > j, 

q~ l \vi)\vj) fti = j. 



A 



Note that when q = 1, we recover the the action defined by 4.4.1. 

In the rest of this thesis we consider the case where V is the representation V 
of U q (d) indexed by the single-box partition A = (1). Using the coproduct struction 
on U q (n), we interpret V® n as a representation of U q (n) as well as H q (n). In order 
for V® n to be a representation of the algebra U q (d) (8) H q (n), their respective actions 
must commute. This is proved by Jimbo in [7], as well as the following result which 
is known as quantum Schur-Weyl duality. 

Theorem 4.4.1. The space V rm as a representation of U q (d) ® H q {n) decomposes 
into irreps in the following formula 

yrm ^ V x (g)R x (4.4.2) 

\\-n,e(\)<d 

Given that the basis of V is indexed by SSYT with a single box, we think of basis 
elements of V as letters, and thus the natural basis vectors of V® n as words of length 
n. This basis can be called the the computational or word basis for V® n . 

With respect to the word basis for V® n and the GTJ and YYH bases for V x 
and R x , respectively, an algorithm carrying out the isomorphism given in 4.4.1 is 
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known as a Schur-Weyl transform. In the next chapter we describe a transform 
which decomposes V® n and prove that it is in fact a Schur-Weyl transform. The 
transform we define has a nice recursive structure, and we prove it has efficient time 
complexity in Chapter 6. 
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Chapter 5 

The Pieri and Schur-Weyl 
transforms 

5.1 Introduction 

In Chapter 4 we presented the representation theories of the quantum group U q (d) and 
the Hecke algebra H q (n), and their correspondence via Schur-Weyl duality. In this 
chapter we present a transform which we prove computes a Schur-Weyl transform. 
In Section 5.2 we define a Pieri transform via only the representation theory of the 
quantum group U q (d). In Section 5.3 we compose Pieri transforms and prove that 
this computes a Schur-Weyl transform. This result is stated without proof in [1], and 
we were unable to find a proof in the literature. Finally in Section 5.4 we look at 
the Pieri and Schur-Weyl transforms in their crystal limits, which ties representation 
theory together with the insertion algorithms seen in Section 3.3. 
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5.2 The Pieri transform 
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A Pieri rule is a formula for decomposing representations of U q (d) <g> U q (d) that take 
the form V x (g) . In this chapter we will only need a Pieri rule for the case when 
m = 1 so that we are decomposing V x ® V. 

Although the inclusion U q (d) C U q (d) (g> is not in general a Gelfand pair, in 
our case of interest the branching rule is multiplicity-free and is given in the following 
theorem. 

Theorem 5.2.1. Given the representation V x and the representation V = of 
U q (d), their tensor product decomposes into irreps according to the following formula. 

V X ®V^ V (5.2.1) 

fi covers X 
£(\)<d 

An algorithm carrying out the isomorphism in equation 5.2.1 with respect to the 
GTJ bases for all U q (d) representations is called a Pieri transform. Interpreting the 
Pieri transform at the level of tableaux, q-insertion becomes relevant. A basis element 
of V x ® V is indexed by an SSYT t and a letter i. The correct way of thinking of 
the Pieri transform is that it q-inserts i into t in superposition, resulting in elements 
which index vectors in for /j, covering A. A matrix entry of the Pieri transform is 
written (s \ t, i) and is non-zero only in the case when s is a result of q-inserting i 
into t. (This follows directly from the behavior of the generators of U q (d).) The Pieri 
coefficients are also called Wigner coefficients. 

In the rest of this section we describe formulas for (s \ t,i). The formulas can 
be found in a variety of sources, in particular in [3]. Similar to the GTJ coefficients 
described in Theorems 4.2.2 and 4.2.3, they are complicated-looking but in fact arise 
from simple combinatorial properties (e.g. axial distances) of the tableaux associated 
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to the basis elements. 

As mentioned above, when computing (s \ t,i), we visualize i as being q-inserted 
in t, activating a sequence of letter bumps. The first letter considered is i, and it 
must bump a larger letter until the last letter bumped is some letter < d. WLOG 
we may assume that i& = d, because if not we may work over the smaller algebra 
U q (ik)- Thus we can visualize a sequence of letters 



influenced by the bumping process. The Wigner coefficients factor into a product of 
reduced Wigner coefficients, one for each letter in 5.2.2. The letter i gets a special 
type of reduced Wigner coefficient, which we'll refer to as type zero. 

Theorem 5.2.2. Suppose i is q-inserted into the SSYT t into the box b in row r. 

The corresponding type zero reduced Wigner coefficient is given by 



Any other letter k in the chain 5.2.2 is assigned a reduced Wigner coefficient of 
type one. Here we have two parameters: the box k inhabits in t in row r\ and the box 
k gets bumped into in tableau s, in row r 2 . Then the type one Wigner coefficients 
are given in the following theorem. 

Theorem 5.2.3. Suppose i is q-inserted into the SSYT t and k gets bumped from 
row ri to r 2 as a result of the insertion. The corresponding type one reduced Wigner 



ii < i 2 < ■ ■ ■ < ik = d 



(5.2.2) 




(5.2.3) 



where are the axial distances with respect to the horizontal strip . 
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coefficient is given by 



W 1 (k; AW) = sgn{n - r 2 )q^-^'\ \ [] + ^ ± 1] II 



[ a j>2 Aj ] 



j<k 



[%> 2 + 1] ;<V-i [a iri + A^ } - Af 1 



(5.2.4) 



where djf are the axial distances with respect to the SSYT t^ and sgn(0) 



Note that the exponent of q given by a rir2 — A^ in theorem 5.2.3 is simply the 
distance from the old box k inhabited in t to the new box k inhabits in s. 

Theorem 5.2.4. The value of the Wigner coefficient (s \ t, i) with associated sequence 
(5.2.2) and with horizontal strips = sh(t^) \ s/i(^' t_1 ^) is given by 



[ s \t^) = w (i)Y[w 1 (i 1 ) 

i=2 



(5.2.5) 



Example 5.2.5. Let t 



so that i = 2 is the letter q- 



and s 

inserted into t to make s. There are two reduced Wigner coefficients, one associated 
to 2, and one associated to the 3 that 2 bumps. 

1 I 1 I o 

J. The 



f 


f 


2 


2 





In the first step, 2 is added onto the second row of the tableau 
residue of this box is zero, X2{t) = 2, and the relevant axial distance is ayi = 3. Thus 
the type zero reduced Wigner coefficient associated to the 2 is given by 



W (2) = q 



-1/2 



In the second step, the 3 is bumped from its original position in the second row 
to its new position in the third row. The axial distance from the old box to the new 
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box is given by 1. The relevant axial distances are an = 2, a 2 3 = 2 and ai 3 = 4. 




Therefore, the Wigner coefficient (s \ t, i) is given by the product 



W (2)W 1 (S) 



M [2] 
[3]v^J 



Theorem 5.2.4 proves that the matrix entries of the Pieri transform decompose 
into a product of reduced Wigner coefficients. Another way of interpreting formula 
5.2.5 is the following recursive version. 



We then define the reduced Wigner transform to be an algorithm computing the d X d 
matrix of reduced Wigner coefficients where t is fixed but % and the nonzero row in 
sh(s) \ sh(t) both vary. 

5.3 The Schur-Weyl transform 

Recall we use the notation V for the d- dimensional irrep of U q (d) indexed by the 
single-box partition A = (1). Thus we think of basis elements of V as just letters and 
basis elements of V® n as words of length n. 

For n = 2, we have V® 2 which can be decomposed with a single Pieri transform 
seen in Section 5.2. For larger n, the Pieri transforms can be composed (or, cascaded) 
to create a transform with input space V® n . For example, when n = 3, we realize a 



(s | t,i) 



Wtid)^- 1 ^ | t^J) i ^ d 



(5.2.6) 



= < 



W (d) 



i = d 
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decomposition of V® n via two sequential applications of the Pieri transform. 
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V® 3 = (V®V)®V 

* (yt 1 ' 1 '®^)®^ 
(v (1,1) ® v) © (V (2) ® y) 

^ (1/(2,1) ^d,l,D) (y(3) 0y (2,l)) 
^(1,1,1) (^(2,D g, y(2,l)) y(3) 

The above description can be read as two sequential q-insertions mapping a word of 
length three to a tableau with three boxes. Note that there are two copies of V^ 2 ' 1 ^ 
in the decomposition, determined by whether the box in the second row was added 
in the first or second instance of the Pieri transform, producing a multiplicity space 
for the irrep V^ 2,l \ As indicated by the notation, we know the multiplicity 
space RC 2 ' 1 ^ is isomorphic to the irrep R^ 2 ' 1 ^ of the Hecke algebra by the Schur-Weyl 
duality theorem 4.4.1. 

By induction, a decomposition of V® n can be achieved by cascading (n — 1) Pieri 
transforms, which we refer to also as a Pieri transform. This transform results in 
a sum of U q (d) irreps V x , whose multiplicity spaces are isomorphic to H q {n) irreps 
R x via the Schur-Weyl duality theorem. However, a priori it is unclear whether the 
change-of-basis achieved by our Pieri transform is in fact identical to the change-of- 
basis required by a Schur-Weyl transform. In the rest of this section we prove that 
cascaded Pieri transforms compute the Schur-Weyl transform up to sign. 

We refer to the GTT basis achieved by (n—1) Pieri transforms as the Pieri basis 
and the basis for Schur-Weyl duality as the Schur basis, as in Section 4.4. 

Theorem 5.3.1. The Schur transform for decomposing V® n by (n — 1) cascaded Pieri 
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transforms is a Schur-Weyl transform as defined in Section 4-4 U P to signs. 

Proof. We first argue the q — 1 case, so that we can work with groups. Beginning 
with the group S{n) xU(d), in both the Pieri and Schur-Weyl decompositions we have 
subgroup flags that reach first U(d) to produce a summand of irreps V x . From this 
point we continue with the subgroup flag defined by the branching rule in theorem 
4.2.1, to yield the standard GTT basis of V x . 

In the Schur basis the subgroup flag is defined by 

U{d) C S{2) x U{d) C S(3) x U{d) C • • • C S{n) x U{d) 

and in the Pieri basis the subgroup flag is defined by 

U(d) C U{d) 2 C [7(d) 3 C • • ■ C U{d) n 

In the Pieri flag, we mean more precisely that within the group U(d) k , the first factor 
of U (d) should act on the first n — k + 1 tensor factors of V diagonally, while the jth 
factor of U(d) for j > 2 should act on the (n — k + j')th tensor factor of V. In other 
words, we can define the desired embedding U(d) k C [7(<i)' c+1 by the map 

(A x id^ 1 ) : U(d) k -> 

where 

A : t/(d) -)• U{df 

is the standard diagonal embedding, and id is the identity. 

We claim that the direct sum decompositions induced by the two subgroup flags 
become equal when they reach U(d) and irreps V x of this group. Since the subgroup 
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flags coincide below U(d), the decomposition must remain equal afterwards. To 
prove the claim, we consider the partially ordered set P of groups S(j) x U(d) k , with 
j + k < n + 1, as shown in Figure 5.3.1. 

P has a unique minimal element, U(d), and n maximal elements. It also has 2 n 
maximal chains that connect U (d) to some maximal element, all of length n+ 1; each 
step of such a chain, in the figure, can either be down or to the right. We claim that 
each of these maximal chains are all locally multiplicity free on V® n , and that the 
partial decompositions all coincide when they reach U(d). 

To see that each chain c C P produces a GTT line basis, we first decompose V® n 
as a representation of a maximal group S(k) x U(d) n ~ k+1 . By 4.4.1 we obtain the 
following isomorphism: 

= [ R x ® V x ) ® V® n ~ k 

\Xhk ) 

This is multiplicity free. Each subsequent step of c is one of the inclusions 

S(j - 1) x U{d) k C Sf(j) x 
S(j) x ^(rf)^ 1 C Sf(j) x U(d) k 

Both of these inclusions are Gelfand pairs by Theorems 4.2.1 and 4.3.1, and the 
structure of irreps of the direct product of two groups. 

To see that the decompositions coincide, we consider two types of moves on chains 
in P: a triangle move that changes the last step between horizontal and vertical, and 
a square move that switches a horizontal and vertical step lower in the chain. The 
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triangle move relates two chains ci, c-i C P that agree except at the three groups: 

S(k) x U{d) n - k C S(fc + 1) x C/(d) n "* 

in 

£>(fc) x f/(rf) n - fc+1 

We claim that the decomposition of V® n ~ k is multiplicity free using the chain C3 = 
c\ Hc2, which begins directly with x U(d) n ~ k . By 4.4.1 and either Theorem 4.2.1 
or Theorem 4.3.1, we obtain 



m—k—l 



/xcoversA 

which is multiplicity free. Since c\ and C2 each yield the same decomposition as C3, 
they yield the same decomposition as each other. 

Likewise suppose that ci,C2 C P differ by a square move: 

S(j) x U{d) k C S(j + 1) x f/(d) fc 

in in 
5(j + 1) x t/(d) fc C S(j + 1) x t/(d) fc+1 

We claim that C3 = ci D C2 is again locally multiplicity free, which implies that c\ 
and C2 must each yield the same decomposition as C3. At the lower left corner, a 
single summand which is an irrep of S(j + 1) x U(d) k+1 will in general have the form 
(R\ <g> Vfj) ® V® fe . Then its restriction to S(j) x U (d) k is multiplicity free, by applying 
Theorem 4.3.1 to R\ and Theorem 4.2.1 to <g> V. 

It is easy to see that all maximal chains in P are connected by square and triangle 
moves. This yields the result when q = 1. 
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[/(d) C 5(2) x [/(d) C 5(3) x U(d) C 

in in in 

[/(d) 2 C 5(2) x U{df C 5(3) x (7(d) 2 C 

in in in 

(7(d) :i C 5(2) x [/(d) :i C 5(3) x U{df C 

in in in 

in in in 

[/(d)'- 2 C 5(2) x U(d) n -' 2 C 5(3) x (7(d)"~ 2 

in in 

[/(d)"- 1 C 5(2) x (7(d)"- 1 

in 

[/(d)" 

Figure 5.3.1: The poset P of group inclusions 

When q is positive (or more generally, when q is not a root of unity) we can follow 
the same argument, except that we replace S(j) x U(d) k by H q (j) <g> U q (d)® k . The 
replacement yields well-defined algebra actions by 4.4.1, and the argument still works 
because it relies on the same multiplicity free structures. 

The line basis agreement extends to unique vector bases up to sign. Because 
the representations considered are all bar representations, the same Gelfand-Tsetlin 
constructions yield unique real line bases. Then, because the representations are 
all ^representations, the real line bases can be refined to bases of real unit vectors. 
These vectors are then unique up to sign. 

□ 



5.4 Pieri and Schur-Weyl in the crystal limit 

In quantum algebra the limit q = is referred to as the crystal limit. The quan- 
tum groups and Hecke algebras we have considered do not have well-defined algebra 
structures at q — 0. However, it can be useful to look at the transforms in the crys- 
tal limit, as they are still linear maps between vector spaces, if not representation 



C S(n - 2) x [/(d) C S(n - 1) x [/(d) C 5(n) x 

in in 
C S{n - 2) x [/(d) 2 C S{n - 1) x [/(d) 2 
in 

c S{n - 2) x [/(d) 3 
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isomorphisms. In this section we describe the behavior of the Wigner coefficients 
defined in the Pieri transform in Section 5.2 in the crystal limit. These theorems can 
be found in paper [3]. 

Theorem 5.4.1. The type zero Wigner coefficient W${r\ A^) is zero in the limit 
q = except when for every j between r and i — 1 we have = X^ ^ . The type 
zero Wigner coefficient Wo(r; AW) is zero in the limit q = oo except when r — 1, and 
in this case it's one. 

Theorem 5.4.2. The type one Wigner coefficient Wi(ri;r 2 ; A^) is zero in the limit 
q = for every choice except when n = r 2 or when n > r 2 and aS^ = ^ for j 
between r 2 and r\ — 1, and in both these cases it's 1. The type one Wigner coefficient 
W / i( r i; r 2; AW) is zero in the limit q = oo for every choice except when r 2 = T\ + 1, 
and in this case it's —1. 

At the level of tableaux, theorems 5.4.1 and 5.4.2 can be interepreted in the 
following way. 

In the limit q = 0, the only time a letter can be inserted into a tableau without 
bumping another letter is when it is added to the very first column that does not 
contain an i. Also, the only time a letter can be bumped, it gets inserted into the 
very next column. This is a description of the dual RSK insertion algorithm. 

In the limit q = oo the only time a letter can be inserted into a tableau without 
bumping another letter is when it can be added to the first row, and the only time a 
letter can be bumped, it gets inserted into the following row. This is a description of 
the RSK insertion algorithm together with a bumping sign described in section 3.3. 

Therefore, in the crystal limits q = and q = oo the transform defined in 5.2.1 
is the dual RSK and RSK insertion algorithm with bumping sign, respectively. In 
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this sense, RSK and dual RSK are 'classical' versions of quantum insertion (not to 
be confused with the other notion of 'classical' in this setting, i.e. q = 1). 

Extending Theorems 5.4.1 and 5.4.2 to cascaded Pieri transforms, we conclude 
that in the crystal limits, the Schur-Weyl transform computes the well-known RSK 
bijections on finite sets described in Section 3.3. 
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Chapter 6 

A quantum algorithm for the 
quantum Schur transform 

6.1 Introduction 

In this chapter we present our main theorem, which is a quantum algorithm for 
computing Schur- Weyl duality. A quantum computer uses basic units of information 
called qubits rather than bits. Since the algebra of qubit spaces is quite different 
than that of bits, the possible types of qubit transformations, and thus the maps 
considered by quantum computers in their calculations, are different as well. This 
in turn influences the types of things that can be calculated within a certain time 
complexity. 

In Section 6.2 we describe the basics of quantum probability, which forms the 
algebraic base for quantum computation. In Section 6.3 we review the basic theory 
of quantum algorithms, and in Section 6.4 we present our main theorem. The methods 
we use in the proof of our main theorem are modeled on those found in the paper 
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[1], where they prove the existence of a Schur-Weyl algorithm for the case q — 1. As 
far as we know, this thesis contains the first instance of an algorithm designed for a 
quantum computer for the purpose of decomposing quantum algebra representations. 



6.2 Quantum probability 

In this section we describe some elements of quantum probability and their relation 
to quantum algorithms. A thorough treatment of this material can be found in the 
book [13]. 

We fix the computational basis of C d to consist of the orthonormal vectors |1), . . . , \d), 
and write all matrices in this basis. The algebra Md is the set of all d x d matri- 
ces in this basis. (The more standard numbering in computer science is from to 
d— 1; we use the mathematicians' numbering which is more standard in combinatorial 
representation theory.) 

When d = 2, M 2 is called a qubit, and represents the probability space corre- 
sponding to a quantum particle with two basis states, such as an electron which can 
be measured as either spin-up or spin-down. For larger d, Md is called a audit. Joint 
systems are constructed by tensoring qudits. So, for example, a pair of qubits is 
M 2 ®M 2 = M 4 . 

The state of a qudit Md is defined to be a positive and normalized dual vector 
on Md- There is an isomorphism between Md and its dual space so we can view dual 
vectors of Md as elements of Md itself. Under this isomorphism, p 6 Md acts as a 
dual vector on Md according to the formula 



Pip) = Tr(pa) 



(6.2.1) 
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The pure states of Md are of the form pu\ = \ifj)(ip\, where \ifj) is a normalized 
vector in C d . Following Equation 6.2.1, its action on an element b G is given by 
the inner product 



Thus pure states are indexed by normalized vectors in C . 

The self-adjoint elements in Md, i.e. those that satisfy a* = a, are called observ- 
ables or measurables. The idempotent observables, sometimes called events, satisfy 
a 2 = a, and are interpreted as measuring whether the qudit is in the state a. By 
the Spectral theorem, any observable a can be decomposed into a sum of idempotent 
observables: 



where a(a) is the spectrum of a, and a\ is the projection operator onto the eigenspace 
defined by A. The Spectral theorem guarantees a choice of projection operators a\ 
which are pairwise orthogonal. The probability that the observable a measures A is 



p m (b)=Tr(\i;)(4)\b) = (4)\b\4)). 



a = 




(6.2.2) 




(6.2.3) 



and the state of the qudit passes to the conditional state p\ = p^, which is defined 



by 




(6.2.4) 



We will mainly use the special case a = J2k \ k){k\, which is a complete measure- 
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ment in the computational basis. Expanding the vector \ip) = Ylk^klk), we have 



3,k 

where pj^ = ip*ipk- Therefore, following Equation 6.2.3, the probability that the state 
of the qudit is measured to be \k) is pk,k = li^kl 2 , and following Equation 6.2.4, the 
state then passes to pk = \k){k\. So, after a complete measurement the state of a 
qudit is defined by one of the basis vectors \k), which can be used as an 'answer' to 
a computational question. 

The state of a qudit can also undergo reversible unitary evolution E given by 
conjugations E(p) = UpU" 1 where U G PSU(d), the space of projective unitary 
maps. If p\^) is a pure state, then 

E (pm) = Up^u -1 = pu\^, 

is also pure and defined by the action of an element in PSU(d). 

Pure states, measurement operators, and unitary evolution are the basic elements 
necessary to define quantum algorithms, which we describe in Secction 6.3. 

6.3 Quantum algorithms 

In this section we review basic quantum algorithms. For a more complete introduc- 
tion, refer again to [13]. 

In the context of quantum computing, unitary operators acting on qudit spaces 
are called quantum gates. In the usual interpretation, a quantum gate acting on m 
qudits can act on n > m qudits by acting by U on m of the qudits and the identity 
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on the other n — m qudits. A quantum circuit is a composition of quantum gates. 
The time complexity of a quantum circuit is the number of quantum gates in its 
decomposition. We require that a family of quantum circuits be uniform, meaning 
there is a classical algorithm to compute the decomposition of the quantum circuits 
into quantum gates. The time complexity of this decomposition algorithm is counted 
toward the total time complexity of the quantum algorithm. 

Because the set of quantum gates is infinite, but a quantum computer would have 
access to a finite number of quantum gates, it's generally not possible to construct 
an exact quantum circuit for calculating a given unitary transformation. In other 
words, some approximation will usually be necessary. 

A finite set of gates G that acts on at most m qudits of size d is called universal 
if it generates a dense subgroup of U(d m ) for some m > 2; it consequently densely 
generates U(d n ) for any n > m. In other words, a set G of quantum gates is universal 
if every unitary operator A G U(d n ) can be approximated by a quantum circuit 
composed of elements from G. 

The Solovay-Kitaev theorem [4] establishes that any operator A E U(d) can 
be approximated by words in a universal gate set with time and gate complexity 
poly(<i, log e _1 ), where e is the error of the approximation. This is provided that the 
matrix entries of the gates and the matrix entries of A can be approximated with 
the same time complexity. Thus, up to polylogarithmic overhead, any universal gate 
set is equivalent to all unitary operators in U(d) whose matrix entries can be com- 
puted quickly. Often this theorem is used for fixed values of d, but the algorithm is 
constructive and it is easy to establish that the gate complexity (and classical time 
complexity to choose the gates) is polynomial in d as well. The Solovay-Kitaev theo- 
rem establishes that this time complexity is independent of the universal gate set up 
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to a polylogarithmic factor (assuming that the matrices of the gates can be computed 
quickly) . 

Both classically and quantumly, an efficient algorithm is one which has polynomial 
time complexity. The class of decision problems with efficient classical probabilistic 
algorithms is called BPP, for bounded-error polynomial time. The class of decision 
problems with efficient quantum algorithms is analogously called BQP, for bounded- 
error quantum polynomial time. The class BQP contains the class BPP, meaning 
that quantum algorithms can efficiently solve any decision problems that classical 
algorithms can solve. Whether or not BQP is strictly larger than BPP remains an 
open problem. However, there are decision problems provably in BQP that are not 
provably in BPP, such as the factoring problem investigated in [16]. 



6.4 Time complexity of Schur-Weyl 

In this section we describe the time complexity of computing first a single Pieri 
transform, and and then the cascade of Pieri transforms which calculates a Schur- 
Weyl transform. 

In Equation 5.2.5, we calculated the matrix entry (s \ t, i) of the Pieri transform 
in terms of reduced Wigner coefficients. Using Formulas 5.2.2 and 5.2.3 with 5.2.5 
allow us to calculate all the matrix entries of the Pieri transform, but at this point 
a direct application of the Solovay-Kitaev theorem is not quite sufficient, since the 
entire Pieri operator has a matrix of size d|SSYT(A)|, which is larger than poly(n, d). 

This is why we use the recursive version of 5.2.5 given by 5.2.6 to define the reduced 
Wigner transform as the dx d matrix of reduced Wigner coefficients where t is fixed. 
As described in [1], given an input t, the coefficients in the reduced Wigner transform 
are calculated by a classical algorithm and then the reduced Wigner transform is 
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calculated by a quantum algorithm. 

Theorem 6.4.1. There is a quantum algorithm for computing the controlled reduced 
Wigner transform with time complexity poly(n, d, log e" 1 ), where e is the desired ac- 
curacy. 

Combining the recursive formula 5.2.6 with Theorem 6.4.1, we obtain the following 
result. 

Theorem 6.4.2. There is a quantum algorithm to compute the quantum Pieri trans- 
form with accuracy e with time complexity poly(n,d, log e _1 ), where e is the desired 
acuracy. 

As proved in section 5.3, a Schur-Weyl transform can be computed up to sign 
using (n — 1) cascaded Pieri transforms. We use this fact combined with theorem 
6.4.2 to calculate the total time complexity of the Schur-Weyl transform given in the 
following theorem. 

Theorem 6.4.3. There is a quantum algorithm for computing the cascaded Pieri 
transform which is a Schur- Weyl transform on V® n with time complexity poly (n, d, log(e 
where e is the desired accuracy. 

Since Theorems 6.4.2 and 6.4.3 are entirely based on unitary operators for any 
real number q > 0, and since unitary groups are compact, we can expect the Schur 
algorithm to have well-defined limits at q = and q = oo, provided that we keep 
q real and positive. This expectation turns out to be correct. As proved in [3] and 
discussed in Section 5.4, in the limit q = the transforms converge to permutations 
matrices, while in the limit q = oo, they converge to signed permutation matrices. 

Therefore, at q = and q = oo our algorithm is a unitary version of the dual RSK 
and RSK algorithm combined with a bumping sign algorithm. 
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Chapter 7 
Conclusion 

We conclude this thesis by discussing some directions for further research into the 
topic of Schur transforms. 

First of all, it is still open whether there is an efficient algorithm for the Schur- 
Weyl transform for any value of q, including q = 1, which is jointly polynomial in n 
and logd, in other words polynomial in the input qubit length n(\ogd). 

A second line of investigation involves finding algorithms for other transforms of 
quantum algebras. In particular, the algorithm in [2] which is an efficient computa- 
tion for the quantum Fourier transform for the symmetric group S(n) could possibly 
be generalized to a quantum Fourier-like transform decomposing the regular repre- 
sentation of the Hecke algebra H q (n). 

The quantum Fourier transform decomposes a representation of S(n) x S(n), while 
the Schur transform decomposes a representation of S(n) x U(d). There is a known 
algebra isomorphism for decomposing an analogous representation of U(d) x U(d) 
known as Howe duality, which is extended to the quantum algebra U q (d) in [18]. 
Whether there is an efficient quantum transform for computing a version of Howe 
duality for any value of q, including q — 1, is open. A possible goal is to generalize 
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all three algorithms (Fourier, Schur, and Howe) into a single algorithm. 

Finally, one may investigate the types of applications g-deformed quantum algo- 
rithms can solve. One of the main reasons that quantum Fourier transforms have 
been studied so extensively is their link to hidden subgroup problems which are in 
turn linked to interesting computational problems such as factoring and graph iso- 
morphism. In [1] some applications of the Schur transform for q = 1 are proposed. 
The types of computational problems quantum algorithms for decomposing quantum 
algebra representations help to solve remain to be investigated. 
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